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ABSTRACT 

We investigate the galactic axes orientations within 18 selected clusters, sub-structures 
of the Local Supercluster. For every cluster we map the parameter An (Flin & 
Godlowski 1986) describing the galactic axes alignment with respect to a chosen clus- 
ter pole, divided by its formal error u(An) (s = Au/a(Au)). The cluster pole co- 
ordinates change along the entire celestial sphere. The resulting maps are analyzed 
for correlations of its maxima with directions from the cluster centre to 1.) the de- 
rived 'physical' cluster poles, 2.) the Local Supercluster centre, 3.) the Virgo A centre 
and 4.) the Earth, i.e. along the line of sight (= LOS). The strong maxima - with 
one exception - exist only for non-spiral (NS) sub-samples, with the maximum well 
correlated with the LOS direction. Another of the studied directions may occur close 
to the maximum only if they are close to the LOS and they do not correlate with 
other features visible on the maps. For clusters with a clearly defined maximum of s 
below 3.0 the conclusion generally does not change. For the spiral (S) sub-samples the 
maps are usually at the random noise level. In these cases a weaker, but still existent 
correlation with the LOS is observed and no other evident correlations are noted. We 
conclude that the strong systematic effect, generated by the process of galactic axis 
de-projection from its optical image, is present in the catalogue data. It can mask 
any weak alignment existing in the analysed clusters. With the use of a simple model 
for the systematic effect we are able to reproduce the main characteristic features of 
the maps for non-spiral galaxies. We note, however, a few clusters showing significant 
differences with respect to this model. 
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1 INTRODUCTION 

Some scenarios of cosmological structures evolution predict 
that the orientation of galactic axes in clusters should prefer 
a certain direction while in the other scenarios galaxies are 
expected to be randomly oriented (cf. Shandarin 1974, Wes- 
son 1982, Silk & Efstathiou 1983, Dekel 1985). Attempts to 
reveal any deviation from isotropy were performed over the 
past years with different, and sometimes contradictory re- 
sults. The work done till 1985 is described by MacGillivray 
& Dodd (1985), who point out that most studies in that 
matter agree with the random galactic distribution within 
the Local Supercluster (LSC) plane, however they suggest 
that galactic planes can be oriented preferentially parallel 
to the LSC plane at some distance from it. In an impor- 
tant paper Kapranidis & Sullivan (1983) analysed samples 
of bright spirals belonging to the LSC and found no strong 
evidence for alignment of these galaxies. 

However, Jaaniste & Saar (1977, 1978) claimed the ex- 
istence of the mean perpendicularity of galactic planes with 
respect to the LSC plane. Since the earlier approaches were 
based mostly on analyses of highly inclined and edge-on 



galaxies, Jaaniste & Saar included all galaxies into consid- 
eration, also the face-on ones. This approach was critically 
discussed and modified by Flin & Godlowski (1986). These 
authors and later Godlowski (1993, 1994) analysed large 
galactic samples within the LSC and came to the conclusion 
that there exists a preferential orientation of galactic planes 
perpendicular to the LSC plane and that there is evidence 
for aligning the galactic rotation axes along the directions 
toward the Virgo cluster centre. The comparison of the align- 
ment properties of spiral and non-spiral galaxies reveals the 
physically unexpected conclusion that spirals exhibit weaker 
alignment. It occurs in spite of the fact that the alignment 
should result from the correlated galactic angular momenta 
in the cluster and the spiral galaxies carry substantial an- 
gular momenta allowing for more precise determinations of 
its' spatial orientations to measure such correlations. 

The recent study of Parnovsky et al. (1994) used both 
the analysis of the galactic position angles and the directions 
of galactic axes obtained from de-projections of their images 
in UGC and ESO catalogues. Additionally, the catalogue of 
flat edge-on galaxies compiled by the same authors was used. 
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They detected an anisotropy with the galactic axes distri- 
bution forming a three-axial ellipsoid, showing an excess of 
about 20% in the direction (4 h -6 h , 20° - 40°), and a deficit of 
about 25% in the direction (13 h -15 h , 30°-40°). The results 
are in general agreement with the earlier result of Fliche & 
Souriau (f990) for orientations of the extended HI galactic 
envelopes and the 'cosmic pole' at (5 h 30 m , 7°) derived from 
distant quasars. Later discussion by Flin (1995) points out 
that the anisotropy direction found by these authors is es- 
sentially the same as the result of Flin & Godlowski (1986). 
However, the detected galactic orientation anisotropy is not 
global, but related to the LSC. 

Several aspects of galactic orientations within separate 
clusters were investigated in a few aspects. There seems to 
exist compelling evidence for strong alignment of the CD 
galaxy with its cluster, and the effect is stronger in more 
elongated clusters (Struble 1987, 1990; Mandzhos 1987; van 
Kampen & Rhee 1990; Trevese et al. 1992, Han et al. 1995). 
Muriel & Lambas (1992) showed the existence of system- 
atic effects in orientations of galaxies with respect to their 
neighbours. For spiral galaxies the effect exists with respect 
to the nearest neighbour, while for elliptical galaxies this is 
also true for all neighbours at the distances smaller than 3 
Mpc h -1 . The study of orientations of weaker galaxies in 
clusters (Trevese et al. 1992) does not reveal any signifi- 
cant alignment, confirming similar earlier reports. However, 
the reports of Godlowski (1993, 1994) based on the anal- 
ysis involving de-projection of galactic rotation axes (Flin 
& Godlowski 1986) reveal the possibility of the existence of 
some systematic differences between different substructures 
of the LSC and suggest a possible 'hedgehog- like' orientation 
pattern. 

The work done so far, claiming the detection or non- 
detection of a galactic alignment in galaxy clusters and other 
LSC substructures, with the discrepant alignment character- 
istics, requires further study and clarification. Since the large 
scale correlations between galactic orientations can be de- 
creased by mixing several differently oriented sub-structures, 
the investigation of separate clusters should reveal a much 
clearer picture of the preferred orientation, if any. Therefore, 
the present work was originally intended to investigate the 
galactic axes orientations within 18 selected clusters, sub- 
structures of the Local Supercluster, based on data from the 
Tully's Nearby Galaxy Catalogue (Section 2). For a majority 
of these clusters the number of galaxies is not very large. For 
this reason we have decided to use the statistical parameter 
of the Fourier test, An, describing the preferential orienta- 
tion of galactic axes with respect to the main axis of the clus- 
ter reference frame (Flin & Godlowski 1986; Godlowski 1993, 
1994; see Appendix B). In Section 3, for every investigated 
cluster, we map the value of An divided by its formal error 
cr(An), with the cluster pole co-ordinates changing along 
the entire celestial sphere. The resulting maps are analyzed 
for correlations of the An/cr(An) maxima with directions 
from the cluster centre to 1.) the derived cluster poles, 2.) 
the Local Supercluster centre, 3.) the Virgo A centre and 
4.) the Earth, i.e. along the line of sight (= LOS). For any 
cluster we divide the full sample according to morphological 
type into spiral (S) and non-spiral (NS) sub-samples. The 
strong maxima of An/cr(Aii), above 3.0, with one excep- 
tion occur only for the NS sub-sample, where the maximum 
is well correlated with the LOS direction. The maxima are 



correlated with other directions when those directions are 
close to the LOS. These directions do not correlate in a clear 
way with other features visible on the maps. For structures 
with a clearly defined Aii/a(An) maximum below 3.0 the 
conclusion generally does not change. For the S sub-samples 
the maps are often at the noise level with the maxima of 
Aii/cr(An) ?s 1.0. In these cases no clear correlation of the 
considered directions and features on maps may be seen. We 
conclude that the strong systematic effect, generated by the 
process of galactic axis de-projection from its optical image, 
is present in the catalogue data. It can mask any existing 
weak alignment in the analysed clusters. The divergences 
from the expected form of the map observed in some cases 
may indicate the existence of non-random galactic distribu- 
tions. The effects however are too weak to be considered 
quantitatively with the present method. In section 4, with 
the use of a simple model for this systematic effect, one 
is able to reproduce the main characteristic features of the 
maps for non-spiral galaxies. A short summary and discus- 
sion is given in the last section (Section 5). The present 
analysis provides a firm argument, that one should be very 
careful when using of methods involving reproductions of 
galactic axes orientations from the shapes of their images in 
statistical investigations. The available data for non-spiral 
galaxies are subject to large systematic errors. The use of 
such data in the alignment analysis for samples distributed 
over large parts of the sky could lead to the detection of 
spurious alignments toward regions containing the highest 
numbers of galaxies. The analysis of galactic position angles 
is free of such systematic effects. 



2 OBSERVATIONAL DATA 

We considered a galactic sample selected from Tully's (1988) 
Nearby Galaxy Catalogue, including LSC galaxies with ra- 
dial velocities - corrected for the solar motion - smaller then 
2800 km/sec. This catalogue provides estimates for galactic 
distances, group membership and the derived galactic axis 
inclination to the LOS. The galactic position angles neces- 
sary in our analysis are not available in Tully's catalogue and 
were therefore taken from Nilson (1973, 1974) and Lauberts 
(1982, 1989). For about 10 percent of the considered objects 
the position angels are not available in any source. However, 
in general they are the nearly face-on galaxies and their po- 
sition angle values have only negligible effects on the derived 
galactic axes orientations. For the numerical computations 
below we have taken these position angles at random from 
the uniform distribution. Tully claims that the data from 
his catalogue is free of the Holmberg effect, thus no further 
correction for that effect is introduced into the analysis. 

It is well known that the galactic brightness is a very 
poor distance indicator for neighbour galaxies. For this rea- 
son Kapranidis & Sullivan (1983) in their study of galac- 
tic alignment obtained large differences between samples se- 
lected according to radial velocities and brightness (see also 
Flin & Godlowski 1986). Therefore knowledge of the group 
membership derived from radial velocities is crucial for a se- 
rious analysis. Under present considerations we have selected 
such defined clusters of galaxies from Tully's catalogue. We 
have taken into account all clusters which consist of at least 
40 galaxies. We analysed separately the samples including all 
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cluster galaxies and two sub-samples selected according to 
the galactic morphological type. We extracted spiral galax- 
ies as one sub-sample (S) and non-spiral, mostly elliptical 
ones as the second sub-sample (NS). The considered clus- 
ters are labeled with Tully's numbers, while the names of 
the widely known ones are given in parentheses: 11 (Virgo), 
12 (Ursa Major), 13, 14 (Coma), 15, 17, 21 (Leo), 22, 23, 
31 (Antila-Hydra), 41, 42, 44, 51 (Fornax-Eridianus) , 52, 
53 (Dorado), 61 (Telescopium) , 64. An orientation sketch of 
the Local Supercluster of galaxies with indicated positions 
of the discussed galactic clusters - sub-structures of the LSC 
- is presented in Fig. 1. Detailed lists of the galaxies in the 
investigated groups can be obtained upon request from one 
of the authors (WG). 

In studies following the Jaaniste & Saar (1977) ap- 
proach a crucial point is to obtain a correct value for the 
galaxy inclination angle i. A formula enabling derivation of 
that angle from the observed axial ratio q, valid for oblate 
spheroids, is provided by Holmberg (1946) as 

cos 2 i = (q 2 - ql)/{l - ql) . (2.1) 

Tully (1988) in the NGC catalogue used this formula with 
qo = 0.2 for obtaining inclination angles. Only in a statisti- 
cally insignificant number of cases, he obtained inclinations 
from other information, such as the form of rings or spiral 
structure. Then, he increased the obtained values of i by 
3° , in accordance with the empirical recipe given by Aaron- 
sen et al. (1980). Let us note that Kapranidis & Sullivan 
(1983) used such a formula (without the mentioned recipe 
'+3°') only for spirals. Heidmann et al. (1971) noticed that 
the value of qo depends on the galactic morphological type 
and consequently should vary for different types of galax- 
ies: for spirals qo should be less than 0.2, while for elliptical 
galaxies they should be significantly higher. For this reason 
Flin & Godlowski (1986) and later Godlowski (1993, 1994) 
considered the 'true' axial ratio qo to depend on the galaxy 
morphological type. However, as we argue later, it may be 
questioned if such simplified approach may work reasonably 
well. 



3 TESTING THE GALACTIC ORIENTATIONS 
IN CLUSTERS 

3.1 The data analysis method 

Let us shortly summarize the data analysis method of Flin 
& Godlowski (1986; see also Godlowski 1993, 1994) used for 
the investigation of possible galaxy alignments in LSC clus- 
ters. In the analysis one uses the supergalactic co-ordinate 
system (L, B, P) with the basic great circle ('meridian') 
chosen to pass through the LSC centre in the Virgo cluster. 
For any galaxy we consider two parameters: the galactic 
position angle p and the inclination angle i. With the use 
of these angels two orientation angles are determined: 5 - 
an angle between the normal to the galaxy and the LSC 
plane, and r\ - an angle between the projection of this nor- 
mal at the LSC plane and the direction toward the LSC 
centre. In the present analysis the method was applied for 
angles defined with respect to the individual galaxy clus- 
ter main plane instead of the LSC plane. The distributions 



of the two angles S and r\ can be analyzed using statisti- 
cal tests from Flin & Godlowski (1986) shortly summarized 
in Appendix B. Within this method based on de-projection 
of galactic images one obtains two possible vectors normal 
to the galactic plane and both solutions are considered in 
the analysis (Jaaniste & Saar 1977). The method was pre- 
viously used by one of us to study the galactic alignment 
inside the whole LSC. However, when we consider individ- 
ual clusters the number of galaxies involved may be small 
in some cases and not all statistical tests described by Flin 
& Godlowski (1986) will work well (e.g. \ 2 t e st requires the 
expected number of data per bin to equal at least 7; see, 
however, Snedecor & Cochran 1967 and Domanski 1979). 
We base our present work on a test involving the Fourier 
coefficient An (cf. Appendix B). The value of An charac- 
terizes the mean alignment of galactic planes with respect 
to the chosen axis of the reference frame. The positive value 
of An obtained during the analysis of the angle 5 appears 
when the galactic rotation axes tend to be perpendicular to 
the cluster plane, and for the negative value these axes are 
preferentially parallel to that planeQ. 

As described in Appendix A, the galactic cluster main 
plane and the respective axis perpendicular to the cluster 
main plane can be derived by assuming that the galactic 
cluster is a three-axial spheroid with Gaussian density dis- 
tribution along each axis. We establish directions of these 
axes in the considered cluster by fitting the three dimen- 
sional the galactic distribution to galactic positions taken 
from Tully's catalogue. As the actual clusters are often quite 
irregular and subject to inevitable galactic distance errors 
the derived axes are rather formal ones and only sometimes 
contain physical information about the cluster gravitational 
potential structure. In the analysis below, we allow for all 
possible orientations of the cluster's main axis along the ce- 
lestial sphere. Of course, this axis is not related to the previ- 
ously mentioned 'physical' axes. For each particular choice 
of the cluster pole the value of the Fourier parameter An 
for the angle S is derived in order to seek the maximum 
value of |An|. Because we change the 'pole' position along 
the whole sphere, there is no need to repeat the analysis for 
the angle r\. Maps of the obtained An are presented in su- 
pergalactic co-ordinates L and B in Fig. 2. More precisely, 
we plot the value of An divided by its standard deviation 
cr(Aii) (Aii/<t(Aii) = s). It is sufficient to present one 
hemisphere only, the second one is obtained reflecting op- 
positely directed poles. Next to each map the considered 
cluster number, the number of considered galaxies and the 
maximum value of |An/(r(Aii)| in the map are given. For 
S and NS galaxies, the respective symbol is positioned af- 
ter the indicated number of galaxies. For presentation in 
Fig. 2 we selected four clusters (11, 15, 22, 23) with differ- 
ent characteristic features, the full set of maps for 18 clusters 
is provided in Appendix C. For each cluster, maps are pro- 
vided for all galaxies and for both S and NS sub-samples. 
The maps are presented in the order of decreasing maximum 



* If one calculates An for the angle r], the value of An < 
occurs when the projection of the galactic rotation axis at the 
cluster's plane is preferentially oriented perpendicular to the zero 
point direction of rj, and for An < the projection tends to be 
directed along this direction. 
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Figure 1. An orientation sketch of the Local Supercluster of galaxies with the positions of galaxy clusters used in our analysis indicated. 
The supergalactic co-ordinates (X, Y, Z) are given. The cluster numbers are given near symbols indicating galaxy clusters and straight 
lines join these symbols with their projections on the LSC plane. Our position is marked with 'E'. 



value of \s\. On each map we denote the crucial directions, as 

seen from the centre or toward the centre ^] of the considered 
cluster: 1.) the three derived cluster poles (cf. Appendix A), 
2.) the direction to the Local Supercluster centre, 3.) the 
direction of the Virgo cluster centre and 4.) the line of sight 
to the Earth (= LOS). Data from the respective fits for 
all 18 clusters is summarized in Table 1, including the clus- 
ter number according to Tully, numbers of ALL, S and NS 
galaxies in the cluster, position of the Earth as seen from the 
cluster centre (Le, Be) and its distance R in Mpc (assum- 
ing Hq — 75 km/s/Mpc), and the position of the parameter 
s maximum on the map for ALL galaxies, (L max , B max ), 
followed by the corresponding value of that parameter. 

3.2 Analysis Results 

In the majority of the considered clusters we observe two 
'hills' of positive values in opposite directions on the celes- 
tial sphere separated by a circular 'valley' of negative values 
(see Fig. 4 in Appendix C). The positive hill is expected to 
appear on the map if the direction of the galactic cluster 
axis is close to the direction of galactic axes alignment. If 

T We chose one of the symmetric directions seen on the presented 
hemisphere. 



we take the considered axis perpendicular to that direction, 
the alignment is mostly perpendicular to the axis and we 
observe negative values of An. For the alignment parallel 
to a single axis in the cluster the depth of the negative valley 
on the map should be smaller than the height of the hills. 
In the case when the galactic axes align uniformly along the 
plane, one expects to find negative holes in the directions 
perpendicular to that plane, separated by a somewhat less 
pronounced positive 'bank' for the cluster axis directions 
chosen in the considered plane. Thus the structure most of- 
ten observed (see cluster 11 at Fig. 2a, Fig. 4) could be pro- 
duced if the galactic axes tend to align along the line joining 
the considered hills and then the observed alignment effect 
could often be interpreted as statistically significant. Below 
we will argue against such an interpretation and we suggest 
that a systematic effect along the line of sight may be to 
a large extent responsible for the structures visible on the 
maps. One should note, however, that the picture with pos- 
itive hills surrounded by negative valleys is not a strict rule 
for all clusters. In some statistically less-significant cases, 
like cluster 23 in Fig. 2b, one has separate negative minima 
surrounded by a circle of positive values. Sometimes the rel- 
ative values of positive and negative maxima do not fit the 
above described picture well (for example cluster 52 in Ap- 
pendix C with nearly equal amplitudes). For cluster 15, on 
the map for ALL galaxies (Fig. 2a) the observed maximum 



© 0000 RAS, MNRAS 000, 000-000 



Galactic axes orientations in galaxy clusters 5 



Figure 2a. Cluster 11, Cluster 15 
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Figure 2b. Cluster 22, Cluster 23 
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Figure 2. Maps of s = An/cr(Aii) versus the chosen cluster pole supergalactic co-ordinates (L, B). On the maps, the results of the 
real data analysis are shown on the left, while the maps obtained from our systematic effect modelling (Section 4) are given on the right. 
For each cluster the maps are presented for ALL cluster galaxies and for S and NS sub-samples. On the map we indicate the important 
directions, as seen from the centre of the considered cluster: 1.) three cluster poles - see text (full star, square and triangle), 2.) the 
direction to the Local Supercluster centre (open circle), 3.) the direction of the Virgo A cluster centre (open square) and 4.) the line of 
sight from the Earth (asterisk). Near each panel we give the cluster number appended for the sub-samples with the respective symbol S 
or NS, the number of galaxies in the (sub-) sample TV and the maximum value of s on the map. Four examples are presented in Fig. 2: 
a.) Cluster 11, and Cluster 15, b.) Cluster 22 and Cluster 23 . 



alignment direction does not correlate with the LOS and in 
cluster 22 we do not see any significant peaks. 

A significant difference is claimed to exist in alignment 
characteristics of the spiral galaxies and non-spirals ones 
(e.g. Flin & Godlowski 1986). Therefore, in the next step 
of the analysis we consider differences between such sub- 
samples in each cluster (Fig. 2, Fig. 4). For NS galaxies one 
generally observes a positive peak strongly correlated with 
the LOS direction, surrounded by the valley of negative s . 
An exception is cluster No. 52, where a pronounced negative 
'peak' is encountered. For spiral galaxies, a different picture 
is common. In general one does not observe the above 3 
sigma alignment effects (i.e. with s > 3.0). Usually the ob- 
tained picture is consistent with the random distribution 
appended with a very weak systematic correction along the 
LOS. An exception is cluster No. 15, where we observe a 
pronounced positive peak far from the LOS, which could be 
an indicator of real alignment. 

In Fig. 2 and Fig. 4, the evident strong correlation of 
the positive maximum with the LOS direction exists for NS 
galaxies. It is clear that the same correlation on the maps 
for ALL cluster galaxies is to a large extent generated by the 
NS component. In addition, because of the special position 
of the Earth with respect to the investigated clusters the 
LOS and the Virgo cluster direction are often close to each 
other (cf. Fig. 1). In the considered clusters, a coincidence 
between the Virgo direction and alignment maximum occurs 
only if the Virgo direction is close to the LOS, and it never 
happens if it is far from the LOS. In a few cases where the 
Virgo direction is close to this maximum, the LOS is gener- 
ally somewhat closer. This fact suggests that the systematic 
error in determining the galactic axis inclination angle is 
present in Tully's catalogue which was used for this analy- 
sis. The effect - turning galactic 'faces' toward the observer 
to produce the observed 'alignment' - is much more pro- 
nounced for NS galaxies, but can also be detected on the 
maps for the S ones. The last statement can be checked by 
inspecting of all maps presented in Appendix C (Fig. 4) for 
the S sub-samples. In the following section we propose a 
simple model which allows to reproduce the characteristic 
systematic effects seen on the maps. 



4 MODELLING THE CONSIDERED 
SYSTEMATIC EFFECT 

In order to quantitatively evaluate the considered sys- 
tematic effect we have performed a simple modelling of 
the anisotropic galactic distribution in order to reproduce 
Tully's catalogue data. Our null hypothesis is that in fact 
we have an isotropic distribution of galaxy angular momenta 
and any preferred orientation arises due to favouring the 



LOS direction in elaboration of the observational data. This 
fact suggests that the effect should be expressed as a func- 
tion of the galaxy inclination angle cosine /i. Thus, we pro- 
ceeded in the following way. For any considered cluster the 
galactic spatial positions were not modified. We considered 
the same distribution of galactic poles for all clusters, con- 
sisting of an isotropic part and a simple anisotropic correc- 
tion along the LOS, proportional to fi (i.e. we consider only 
the first two terms of the distribution expansion in powers of 
(i, an example with included terms oc fi 2 is analyzed in Ap- 
pendix D). The assumed normalized probability distribution 
has the form 



+ 



2A 



1 + A 1 + A 



(4.1) 



H is the cosine of the angle between the LOS and the galac- 
tic pole and the distribution exhibits rotational symmetry 
along the LOS. For each galaxy in the analysis one uses 
both solutions with the positive and the negative value of 
fx, but, due to symmetry, the range < Li < 1 is considered 
in Eq. 4.1 . The anisotropic part has a simple form with the 
(positive or negative) amplitude parameter A. 

If the effect of anisotropy is systematic it should influ- 
ence all galaxies in the catalogue, not only those contained 
in the discussed clusters. Therefore, for determination of the 
amplitude A we considered all LSC galaxies from the Tully's 
catalogue. ALL, S and NS galaxies were considered sepa- 
rately and the derived amplitude A is appended with the 
respective index ALL, S or NS. Let us note that with thus 
obtained A it was possible to derive the systematic effect for 
the whole LSC to have the similar An as the value derived 
from the original data. For ALL galaxies the best fit gives 



Aal 



0.7 



(4.2) 



An integral amplitude of the anisotropic part is 
f i 



/ 

Jo 



2A AL LfJ-/{l + A A LL)dfi = 0.41 



A fit quality of the data from the Tully catalogue for ALL 
LSC galaxies can be evaluated by an inspection of Fig. 3 for 
both galactic orientation angles 5 and r\. In the figure, we 
presented the catalogue data superimposed over the 'theo- 
retical' distribution derived as a mean from 1000 simulated 
catalogues using the distribution (4.1) with the amplitude 
A value given in (4.2). 

The discussion in the previous section notes a signif- 
icant difference between distributions for spiral and non- 
spiral galaxies. For the spiral galaxies our fitting procedure 
yields a much smaller value of the amplitude A 
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Table 1. A list of analysed clusters: the respective numbers (TV, N$, Ajvs), cluster positions (L c , B c , R), the LOS position at the 
presented map (Le,Be), the maximum s position at the map (L max , B max ), and the value of this maximum 
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Figure 3. A fit of the real data to the theoretical model for 
ALL galaxies in the LSC. The fitted angle is: a.) <5 , b.) r\. For 
each galaxy both possible orientations of the axis ('solutions') are 
included. 



while for the non-spirals 

A NS = 20.0 



(4.4) 



As = 0.15 



(4.3) 



One should remember that for the spiral galaxies some low 
brightness face-on spirals may be missing and the obtained 
As may be underestimated. 

With the use of the derived average values for A (4.2-4), 
in Fig. 2 and in Appendix C we plotted the maps obtained 
for the considered clusters with the model distribution (4.1), 
to be compared with the real data. In each case we derived 
only one model map, and this map - including the statistical 
fluctuations - is presented. One may note that our simple 
model reproduces the original data for clusters with the large 
maxima of An/a(An) relatively well and thus provides a 
strong argument for the presence of the suggested system- 
atic effect. For spiral galaxies the maps are often at the one 
sigma level, but quite often the simulated map is similar to 
the real one with the LOS close to the maximum on the 
map. However, in some cases noticeable deviations between 
the simulated and the real maps exist. In some clusters the 
expected significant maximum along the LOS does not exist 
(cf. Cluster 22 or 15 at Fig. 2). Also the depth of a negative 
valley surrounding the maximum changes in some cases sig- 
nificant in more than expected for the uniform and isotropic 
galactic distribution and observed maximum is much higher 
then expected from the model (cf NS Cluster 11). In our 
opinion, the obtained large value of Ans proves that Tully's 
procedure for derivation of the galactic axis inclination pro- 
vides values not entirely related to the real galactic orienta- 
tions. Probably his use of a mean galactic ellipsoid form for 
such galaxies is in appropriate. 

Finally, we would like to stress that the presented model 
is not expected to reproduce all details seen on the maps. 
Beside the statistical fluctuations, the observed deviations 
between real and model maps can reflect not only the actual 
galactic alignments, but also the fact that the systematic 
error in the catalogue may have a somewhat different form 
than the simple test distribution introduced in the equation 
(4.1) and/or it can arise from fine galactic morphological 
type variations in individual clusters, leading to a slightly 
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different systematic effects. An example of a more thorough 
analysis of the problem is presented in Appendix D for the 
Virgo cluster. 



5 FINAL REMARKS 

We studied the galactic axes alignments in 18 clusters from 
the LSC with the method based on the analysis of de- 
projected galactic axes orientations. The considered statis- 
tical parameter An describes the galactic axes alignment 
along the given direction. Mapping of this parameter in all 
directions on the sky shows that the most evident alignment 
occurs along the line of sight, at least in the majority of sta- 
tistically significant cases. We interpret this finding as the 
result of the systematic error of the galactic axis determina- 
tion in the Tully (1988) catalogue. 

The anisotropy modelling of the previous section shows 
the possibility of obtaining a reasonable fit to the data with 
the simple model (4.1) and the same value of the anisotropy 
parameter A for all clusters. One may note that the intro- 
duced anisotropy is proportional to the galactic axis pro- 
jection on the LOS, i.e. the derived galactic distribution 
underestimates the number of edge-on galaxies. In Tully's 
catalogue this effect is very strong for non-spiral galaxies 
(Ans = 20), but it can also be noticed in the spiral galax- 
ies' sub-sample, where As = 0.15. The analysis of the maps 
for spiral galaxies, where the systematic effects are much 
weaker, do not reveal any meaningful correlation of the con- 
sidered directions - other than the LOS - with the maps' 
characteristic features. Possible problems with de-projection 
of elliptical galaxies was noted by some authors previously 
(cf. Kapranidis & Sullivan 1983). In the present paper we 
provide convincing proof for such an effect to exist in Tully's 
catalogue. In our opinion the discussed LOS effect arises dur- 
ing computation of the inclination angle, where Tully takes 
the 'true' galaxy axial ratio equal to 0.2 for all galaxies, 
independently of their morphological type^]. 

Let us mention that opposed to our finding Bahcall et 
al. (1990) claim that no excess of face-on ellipticals and the 
excess of edge-on spirals exist for the UGC galaxies. The last 
effect arises probably due to the number of missing face-on 
spirals. Our positive fit for As > suggests that either there 
are very few missing face-on spirals in the Tully catalogue, 
or the systematic de-projection errors are significant enough 
for spiral galaxies to compensate for the lack of some face-on 
galaxies (we do not consider the possibility of a real align- 
ment along the elongated structure of the LSC here). 

Let us mention that some previous positive detections 
of galactic axes alignment can be dominated by similar sys- 
tematic effects introduced by the methods involving de- 

t We also performed (Godlowski, in preparation) the analysis of 
galactic orientations in the UGC and the ESO catalogues (Nilson 
1973, Lauberts 1982). In order to derive galactic spatial orien- 
tations, the Hcidmann et al. (1971) recipe for the galaxy shape 
determination and the standard Holmberg formula for the inclina- 
tion angle were used. The model (4.1) was fitted to thus obtained 
data. The obtained values of A were significantly smaller than the 
ones given in (4.2-4). This result confirm our conclusion that the 
de-projection procedure is responsible for the observed systematic 
effect. 



projection of galactic images. In particular, in our analy- 
sis, the LOS effects are much larger than any actual sig- 
nal, yielding a spurious galactic alignment with the direc- 
tion to the cluster. When analyzing larger samples covering 
a substantial part of the sky - with objects distributed in a 
non-uniform way - the effect of this type should introduce 
galactic axes alignment along the direction pointing toward 
the region containing the largest concentration of galaxies. 
In a similar manner, a deficit of alignment would be pro- 
duced in the directions of less populated parts of the sky. 
In our opinion there is no safe approach, excluding the pos- 
sibility of the LOS systematic effects, when attempting to 
derive galactic spatial orientations for large galactic sam- 
ples. The alignment studies involving the galactic position 
angles should not be influenced by any systematic effect of 
the kind discussed in the present paper if the gravitational 
lensing does not play a role. 

Of course, our analysis does not rule out the possibility 
of some real weak galactic alignments within galaxy clus- 
ters. The exceptions from the general picture involving the 
systematic effect mentioned in the text provide some weak 
evidence for the real alignments existing in a few clusters. 
For example the pronounced positive peak for spiral galaxies 
occurring far from the LOS in Cluster 15 may be suspected 
to represent such a real alignment, see also the discussion 
of positive and negative maxima in section 3. However, one 
should be aware of the fact that the statistical significance 
of any such exception may be much lower than the value 
An/cr(Aii) provided at the map due to the large number 
of tested directions, i.e. a large number of statistical trials. 
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APPENDIX A. FITTING OF THE GALAXY 
CLUSTER AXIS 

We use the supergalactic co-ordinate system (L, B, P), 
where P is the supergalactic position angle defined by Flin & 
Godlowski (1986). The co-ordinates of the supergalactic pole 
in the equatorial system 01950 = 285.5 and <5igso = +16 are 
taken from Tammann & Sandage (1976). In our system the 
basic great circle ('meridian') of the supergalactic system is 
chosen in such a way that it passes through the Virgo clus- 
ter centre with co-ordinates 01950 = 186.25, ($1950 = +13.1 
(Tammann & Sandage 1976). In Tully's (1988) catalogue, 
galactic distances and their membership in the correspond- 
ing groups are given. With the use of this data we derive 
positions for cluster's galaxies in an orthogonal reference 
frame. To the obtained distribution we fit the triaxial ellip- 
soid with the Gaussian density distribution along each axis. 
Through the orthogonalization procedure of the covariance 
matrix of a considered distribution, we derive the directions 
of the main axes of the considered cluster. 
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APPENDIX B. THE APPLIED STATISTICAL 
PROCEDURES 

To check the distribution of galactic orientation angles (5, rf) 
we applied statistical tests originally introduced to this prob- 
lem by Hawley & Peebles (1975), later improved by Kindl 
(1987) and described in detail by Godlowski (1993, 1994). 
Below, a short summary is presented of the x^-test, the 
Fourier test and the autocorrelation test. 

Let N denote the total number of solutions for the 
galactic axes (two solutions for any galaxy) in a considered 
cluster, N k - the number of galaxies with orientations within 
the fc-th angular bin, iVo - the mean number of galaxies per 
bin and, finally, No, k - the expected number of galaxies in 
the fc-th bin. In the present derivations we adopted, in most 
cases, a division for n — 36 bins of equal width. As a check, 
in a few cases we repeated the derivations for different values 
of n but no significant differences appeared. 

The presented in section 4 fits of the anisotropy pa- 
rameter A for all LSC galaxies were obtained with the use 
of the statistical tests described below. The x 2 -test of the 
distribution involves the value 



X 



E 



(N k -N , k ) 2 
N ,k 



(51) 



It is a subject of the minimalization in model fitting, with 
the use of the distribution provided by the model, No,k- In 
the present paper this test was used to fitting the parameters 
A from all LSC galaxies and in Appendix D for the analogous 
fit for the Virgo cluster. The main statistical test used in the 
present paper is the Fourier test involving the first Fourier 
mode only. The actual distribution Nk is approximated by 



N k =A?o,fc(l + Aiicos26) fe + A2isin26y 
The coefficients A;i (i = 1, 2) are given as: 



(52) 



P(> A 1 )=exp(-JiV A?) , 
while the standard deviation of this amplitude is 



^ = ( Jb) 



1/2 



(57) 



(58) 



From the value of An one can deduce the direction of the 
departure from isotropy. If An < 0, then, for 6 = S+-k/2, an 
excess of galaxies with rotation axes parallel to the cluster 
plane is present. For An > the rotation axes tend to be 
perpendicular to the cluster plane. 



APPENDIX C. MAPS OF THE PARAMETER S 
FOR 18 CLUSTERS 

The maps of s = Au/u(An) in supergalactic co-ordinates 
(L, E) are presented for all considered clusters (Fig. 4). For 
each cluster, the maps for real data are presented above the 
respective maps obtained from our systematic effect mod- 
elling. The maps are given for ALL cluster galaxies and for 
S and NS sub-samples. On each map we indicated the im- 
portant directions, as seen from the centre of the considered 
cluster: 1.) three cluster poles (full star, square and trian- 
gle), 2.) the direction to the LSC centre (open circle), 3.) the 
direction of the Virgo A cluster centre (open square) and 4.) 
the LOS (asterisk). Near each panel we give the cluster num- 
ber appended for sub-samples with the respective symbol S 
or NS, the number of galaxies in the (sub-) sample N and 
the maximum value of s on the map. 



An 



ELi^ -^0,*) sin20 fc 



ELi^sin^ 
with the standard deviation 



(53) 



(54) 



-1/2 



cr(Aii) = N , k cos 2 26 k 



2 X 1 / 2 



(nNo) 



(55a) 



-1/2 



2 \V2 



a(A2i)=(£iVo, fe sin 2 20 fc j = ( — ) 
The probability of the amplitude 

Ai = (A n +A 2 1 ) 1/2 

being greater than a certain chosen value is given by the 
formula: 



(556) 



(56) 



APPENDIX D. DIFFERENCES BETWEEN THE 
REAL MAP AND THE MODEL 

As an interesting exercise one could apply the discussed 
model to any particular cluster in order to remove the influ- 
ence of the systematic effect. Of course galaxies of different 
morphological types are included in varying proportions in 
individual clusters and therefore the model anisotropy am- 
plitude parameter A may vary among them, deviating from 
the mean value derived by us. Also, the simplified model 
claiming that the systematic effect is simply due to the term 
proportional to the galaxy inclination cosine may be a sub- 
stantial oversimplification in some clusters. Such effects are 
visible by eye among the maps of Fig. 4. However, as even 
a very imperfect removing of the LOS effect will not erase 
the off-axis alignments, we decided to test in one example 
if any remaining structure coincide with other characteris- 
tic directions of the LSC. When dealing with an individual 
cluster we fit an axi-symmetric LOS anisotropy model in- 
volving the next quadratic term of the anisotropy expansion 
in fi. Such a procedure enables a better subtraction of the 
LOS effect, and thus gives us a chance to look for any galaxy 
alignment non- coincident with this axis. Thus, instead of the 
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Figure 4a. Clusters 11, 12, 13, 14 
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Figure 4b. Clusters 15, 17, 21, 22 
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Figure 4c. Clusters 23, 31, 41, 42 
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Figure 4d. Clusters 44, 51, 52, 53 
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Figure 4e. Clusters 61, 64 



distribution (4.1) we use one involving two first terms of the 
expansion of F(fi) in powers of /J^. 

F{fl) ^ l + A + B + 1 + A + B ^ l + A + B tl ■ 

(D.l) 

We fit this model to the data for our cluster with the largest 
number of galaxies, the Virgo cluster (Cluster 11), in the 
following way. For any set of parameters A and B from the 
ranges < A < 2Alsc and — Alsc < B < Alsc we de- 
rived an axi-symmetric model distribution for test galaxies 
placed in the positions of the real ones {Alsc is the ampli- 
tude fitted within the model (4.1) to the respective - ALL, 
S, NS - data derived in section 4 for all LSC galaxies). For 
any such galaxy distribution, a map of the alignment coef- 
ficient An /<r( An) was provided and subtracted from the 
map for the real cluster. The solutions were selected pro- 
viding the smallest residuals on the map (Fig. 5). On the 
left side of the figure we presented the maps for real cluster 
while on the left side we presented the residuals. Of course, 
our maps generated with the use of the model are subject to 
statistical fluctuations and the described procedure provides 
only a good approximation to the best set (A, B). 

§ One should note, that the data for non-spiral galaxies in the 
Virgo cluster does not allow for a fit with the model (4.1) provid- 
ing the LOS effect amplitude as large as the real one. 



Inspection of Fig. 5 shows that subtraction of the model 
distribution does not leave any significant residues for all 
galaxies (A = 0.18, B = 0.37) and for spiral galaxies (A = 0, 
B — 0.13), when a rather weak 2-a alignment for NS galaxies 
(A = 1.6, B = 2.8), non-coincident with any of the tested 
directions was found. In this last case the significance of the 
residual alignment is even smaller than the cited value due 
to the large number of directions tested. Without a much 
more thorough study one is unable to state firmly if any 
real alignment exists in this cluster. 
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